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1 Introduction 

In recent years quaternion matrix equations have been investigated by many 
authors (see, e.g., [1]-[16]). For example, Jiang, Liu, and Wei pQ studied the 
solutions of the general quaternion matrix equation AXB — CYD = E, and 
Liu |3j studied the least squares Hermitian solution of the quaternion matrix 
equation (A ff XA, B^XB) = (CD), Wang, Chang, and Ning [10] derived 
the common solution to six quaternion matrix equations, etc. 

However the understanding of the problem for determinantal representing 
the least squares solutions of the quaternion matrix equations, in particular 

AX = B, (1) 
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XA = B, 
AXB = D 



(2) 
(3) 



has not yet reached a satisfactory level. The reason was the lack of appro- 
priate noncommutative determinant. Many authors had tried to give the 
definitions of the determinants of a quaternion matrix, (see, e.g. [T7]-|23]). 
Unfortunately, by their definitions it is impossible for us to give an determi- 
nant representation of an inverse matrix. 

But in [21], we defined the row and column determinants and the double 
determinant of a square matrix over the quaternion skew field. As applica- 
tions we obtained the determinantal representations of an inverse matrix by 
an analogue of the adjoint matrix in [23] and the Moore-Penrose inverse over 
quaternion skew field EI in [21]. In [16] we also gave Cramer's rule for the 
solution of nonsingular quaternion matrix equations (jTJ), 02]) and within 
the framework of the theory of the column and row determinants. In [131 E], 
the authors obtained the generalized Cramer rules for the unique solution of 
the matrix equation ([3]) in some restricted conditions within the framework 
of the theory of the column and row determinants as well. 

In this paper we aim to obtain explicit representation formulas (analogs 
of Cramer's rule) for the minimum norm least squares solutions of quater- 
nion matrix equations (JT]), 02]) and (j3]) without any restriction. The paper 
is organized as follows. In Section 2, we start with some basic concepts and 
results from the theory of the column and row determinants which are nec- 
essary for the following. The theory of the column and row determinants 
of a quaternionic matrix is considered completely in [23]- In Section 3, we 
give the theorem about determinantal representations of the Moore-Penrose 
inverse over the quaternion skew field derived in [22]. In Section 4, we obtain 
explicit representation formulas for the minimum norm least squares solu- 
tions of quaternion matrix equations (jTJ, ([2]) and ([3]). In Section 5, we show 
a numerical example to illustrate the main result. 

2 Elements of the theory of the column and 
row determinants. 

Throughout the paper, we denote the real number field by M, the set of all 
m x n matrices over the quaternion algebra 

H = {a + a\i + 0,2] + a^k \ i 2 = j 2 = k 2 = — 1, a Q , ai, a 2 , G K} 
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by H mxn and its subset of matrices of rank r by H™ xn . Let M (n,H) be the 
ring ofnxn quaternion matrices. 

The conjugate of a quaternion a = ciq + aiz + a2j + CX3/C G HI is defined by 
a = a Q — a\i — a 2 j — a 3 k. The Hermitian adjoint matrix of A = (a^) G H nxm 
is called the matrix A* = (a* 4 ) if a*, — ~a~n for all % — 1, . . . , n and 

V l / my.n % 3 ■> 

j — 1, . . . , m. The matrix A = (a^) G HP xm is Hermitian if A* = A. 
Suppose 5 n is the symmetric group on the set /„ = {1, . . . , n}. 

Definition 2.1 The ith row determinant of A = (a^) G M(n,M) is defined 
by 

rdetjA = ( — 1) a ii kl a 'i kl i kl + 1 ■ ■ ■ a i kl+h i ■ ■ ■ a i kr i kr+1 ■ ■ ■ a i kr+ i r i kr 

creS n 

for all i = 1, . . . , n. The elements of the permutation a are indices of each 
monomial. The left-ordered cycle notation of the permutation a is written as 
follows, 

a — (i ifcjifcj+i • • • ifci+h) fe^fa+i • • • ik 2 +h) ■ ■ ■ (4 r 4 r +i • • • h r +i r ) ■ 

The index i opens the first cycle from the left and other cycles satisfy the 
following conditions, ik 2 <%<... < ik r and i^ t < ik t + s for all t = 2, . . . , r 
and s — 1, . . . , l t . 

Definition 2.2 The jth column determinant of A = (ay) G M(n, H) is 
defined by 



cdetj A. ^ ^ ( 1 ) a ; . ■ ■ ■ o n . . .( i: j ! t . . . ( 1 i j f ('j 



n—r 

n '3k r 3k r +lr ' ' ' %+1'ir ' ' """JJfci+ii ' ' ' Uj 3k 1 + l3k 1 u '3 kl 3 

r£S n 



for all j = l,...,n. The right-ordered cycle notation of the permutation 
t G S n is written as follows, 

T = {jkr+lr ■ ■ -jkr + likr) ■ ■ ■ (jk 2 +h ■ ■ ■ j k 2 + lj k 2 ) (jfcl+h ■ ■ -jkl + ljklj) ■ 

The index j opens the first cycle from the right and other cycles satisfy the 
following conditions, j k2 < j k:i < . . . < j kr and j kt < j kt+s for all t = 2, . . . , r 
and s — 1, . . . , l t . 

The following lemmas enable us to expand rdetj A by cofactors along the ith 
row and cdetjA along the jth column respectively for all i,j — 1, . . . ,n. 
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Lemma 2.1 Jlffi Let f?jj be the right ij-th cofactor of A E M(n, H) ; i/icrf 

n 

is, rdetj A = a ij ' ^ij / or all i = 1, . . . ,n. Then 

J -rdetj-A^a.i), % ^j, 



rdet fc A", i = j, 



where A" (a j) is obtained from A fry replacing the jth column with the ith 
column, and then by deleting both the ith row and column, k = min {/„ \ {i}}. 

Lemma 2.2 Lei be the left ij-th cofactor of A E M (n, EI), that is, 

n 

cdetj A = Lij ■ a^j for all j = 1, . . . , n. Then 



i=i 



-cdcl, A/;' (a ; .) . i^j, 
cdetfc A J J , % — j, 



where A| J (aj.) is obtained from A by replacing the ith row with the jth row, 
and then by deleting both the jth row and column, k = min{J n \ {j}}. 

Since these matrix functionals do not satisfy the axioms of noncommu- 
tative determinant, then they are called "determinants" shareware. But by 
the following theorems we introduce the concepts of a determinant of a Her- 
mitian matrix and a double determinant which both satisfy the axioms of 
noncommutative determinant and can be expanded by cofactors along an 
arbitrary row or column using the row-column determinants. This enables 
us to obtain determinantal representations of the inverse matrix. 

Theorem 2.1 [2^] If A = (ay) E M(n,H) is Hermitian, then rdetiA = 
■ • • = rdet n A = cdetiA = • • • = cdet„A E R. 

Remark 2.1 Since Theorem \2.1\ we can define the determinant of a Her- 
mitian matrix A E M (n, H) putting for all i = 1, . . . ,n 

det A := rdetj A = cdetj A. 

Theorem 2.2 \2~1$ If for a Hermitian matrix A E M (n,M), 

det A ^ 0, 
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then there exist a unique right inverse matrix (LA) -1 and a unique left in- 
verse matrix (LA) -1 of A, where (LA) -1 = (LA) -1 =: A -1 ; and they 
possess the following determinantal representations 



(4) 







(Ru 


R21 ■ 


" Rnl^ 


(LA) -1 = 


1 


Rl2 


R22 ■ 


Rn2 


det A 












\Rln 


R211 " 


Rnn / 






/L n 


L2I " " 


L n l^ 


(LA)" 1 = 


1 


Ll2 


L22 • ■ 


L n2 


det A 












\Lln 


L2n • • 


L nn J 



(5) 



where Rij, are right and left ijth cof actors of A respectively for all i,j = 
1, ...,n. 

Theorem 2.3 §///AeM(n,H), then detAA* = det A* A. 

According to Theorem 12 .31 we introduce the concept of a double determinant. 



Definition 2.3 The determinant of the corresponding Hermitian matrices, 
(A* A or AA*), of A G M(n, H) is called its double determinant, i.e. ddetA := 
det (A* A) = det (AA*). 

Suppose A lJ denotes the submatrix of A obtained by deleting both the 
ith row and the jth column. Let a.j be the jth column and a«. be the ith row 
of A. Denote by a* and a* the jth column and the ith row of a Hermitian 
adjoint matrix A* as well. Suppose A j (b) denotes the matrix obtained from 
A by replacing its jth column with the column b, and Aj. (b) denotes the 
matrix obtained from A by replacing its ith row with the row b. 

We have the following theorem on the determinantal representation of 
the inverse matrix over H. 

Theorem 2.4 [2J^] The necessary and sufficient condition of invertibility of 
A e M(n, H) is ddetA ^ 0. Then there exists A" 1 = (LA) -1 = (L!A)~\ 

where 

/Ln L21 . . . L n i\ 

,-i -1 1 



(LA)- 



(A*A) _1 A* 



ddetA 



L12 L22 • • • Ln2 
yLi n \^2n ■ ■ ■ ^nnj 



(6) 
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(RA)' 1 = A* (AAT 1 



/Mil M21 

&12 M 2 2 



ddetA* 



-»2 



and 



Ly = odet i (A*A)^(a: 4 ), 
/or a// i, j = 1, n. 



\Mi n M 2ra ... M. 
: rdet^AA*),. (a*) 



Remark 2.2 Since by Theorem\2.4 



ddetA = cdetj (A* A) = V^Ly ■ a y -, ddetA = rdetj (AA*) = V" 



(7) 



for all j = 1, . . . , n, then Ly called the left double ijth cof actor and 
called the right double ijth cof actor for the entry o/AgM (n, H). 



/.s 



3 Deter minantal representation of the Moore- 
Penrose inverse. 

For a quaternion matrix A 6 jj mxn ? a generalized inverse of A is a quaternion 
matrix X with following Penrose conditions 

1) (AA+)* = AA+; 

2) (A+A)* = A+A; 

3) AA+A = A; 1 ] 

4) A+AA+ = A+. 



Definition 3.1 For A e M mxn , X G H nxm zs sazd to fre a gener- 
alized inverse of A if X satisfies Penrose conditions in $3j). We 
denote the X fry A^- 7 '--) and the set of all A^'-^ by A{i,j, ...}. 

By [26] we know that for A e H mxn , the A{«, j, ...} exists and the A^' 2 ' 3 ' 4 ) 
exists uniquely. The matrix A^ 1 ' 2,3 ' 4 ) is called the MoorePenrose inverse of 
A and denote A+ := A^ 1 ' 2 ' 3 ' 4 ). 

In [25] the determinantal representations of the Moore-Penrose inverse 
over the quaternion skew field was derived based on the limit representation. 
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Theorem 3.1 If A G M_ mxn and A + is its Moore-Penrose inverse, then 
A+ = limA* (AA* + aiy 1 = lim (A* A + al)" 1 A*, where a G M + . 

a— >0 a— s-0 

Corollary 3.1 If A £ H mXn ; then the following statements are true. 

i) If rank A = n ; then A+ = (A*A) _1 A* . 

ii) //rank A = m, then A + = A* (AA*) -1 . 
Hi) If rank A = n = m ; then A + = A -1 . 

We shall use the following notations. Let a := {ai, . . . , C {1,..., m} and 
(3 := . . . , /3 fe } C {1, . . . , n} be subsets of the order 1 < /c < min {m, n}. 
By A^ denote the submatrix of A determined by the rows indexed by a and 
the columns indexed by ft. Then A" denotes the principal submatrix deter- 
mined by the rows and columns indexed by a. If A G M (n, H) is Hermitian, 
then by |A^| denote the corresponding principal minor of det A. For 1 < k < 
n, denote by Lfc >n := {a : a = (a±, . . . , o^) , 1 < «i < . . . < ak < n} the col- 
lection of strictly increasing sequences of k integers chosen from {1, . . . , n}. 
For fixed z G a and j G (3, let I r , m {i} '■= { ot : a G £ r , m , z G a}, </ r ,n{j} := 

{P:PeL rin ,jep}. 

The following theorem and remarks introduce the determinantal repre- 
sentations of the Moore-Penrose inverse which we shall use below. 

Theorem 3.2 [25] If A G W? xn , then the Moore-Penrose inverse A + = 
(a^) G H nxm possess the following determinantal representations: 

E cdet^A'A).^)) J 

4 = ^ r ,{i} _ — , (9) 



E 



f A* A) i 



or 



E rdet J ((AA*),.(a*))- 

< = " g7r ' mm ^ , , . . r , . (10) 

13 E l(AA*)«| 



/or aZZ z = 1, n, j = 1, m. 
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Remark 3.1 //rank A = n, then by Corollarv \3.1\ A + = (A* A) 1 A*. Con- 
sidering (A* A) -1 as a left inverse, we get the following representation of A + : 



1 ^cdet 1 (A*A). 1 (a* 1 ) ... cdet 1 (A*A). 1 (a*J 

x cdet n (A*A). n (a* 1 ) ... cdet„(A*A). n (a* m ) 
If m > n, then by Theorem \3.S\ for A + we have |PP as well. 



ddetA 



Remark 3.2 //rank A = m, then by Corollary CO A + = A* (AA*) \ 

Considering (AA*) -1 as a right inverse, we get the following representation 
ofA + : 

1 / rdeti(AA*) x . (a* ) ... rdet m (AA%. (a* ) 

ddetA \rdet!(AA*)i. «) ... rdet m (AA*) m . «. 
If m < n, then by Theorem \3.S\ for A + we also have (Q7J). 



4 Explicit representation formulas for the 
minimum norm least squares solutions of 
some quaternion matrix equations 

Denote by ||A|| the Frobenius norm of the quaternion matrix A £ jj mxn . 

Definition 4.1 Consider a matrix equation 

AX = B, (13) 

where A G H mxn ,B G M. mxs are given, X G H nxs is unknown. Suppose 

H R = {X|X G H nxs , ||AX - B|| = min}, 

Then matrices X G M. nxs such that X G Z/r are called least squares solutions 
of the matrix equation [TB]) . If^Kis — mm xgff fi ||X||, then X^5 is called the 
minimum norm least squares solution of < f73j) . 

If the equation ffl3|) has no precision solutions, then X i5 is its optimal ap- 
proximation. 

The following important theorem is well-known. 
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Theorem 4.1 (JM/) The least squares solutions of /il3\) are 

X = A+B + (I n - A+A)C, 

where C G M nxs is an arbitrary quaternion matrix and the minimum norm 
least squares solution is X^s = A + B. 

We denote A*B =: B = (fey) G H nxs . 

Theorem 4.2 (i) //rank A = r < m < n, then for the minimum norm 
least square least squares solution ~K LS = (xy) G H nxs o/ / f73j) for all 
i = 1 , . . . , n, j = 1 , . . . , s, we have 



E cdeti((A*A) A (b, 

/3eJr,n{«} 



E 



fA*A x 



(14) 



(mJ //rank A = n ; t/ien /or X L5 = (x^) G H nxs of /or a// z = 1, ...,n, 
j = 1, s, we have 



cdet,(A*A) 



X{j 



ddetA 

where h.j is the jth column o/B for all j = 1, s. 



(15) 



Proof, (i) If rank A = r < m < n, then by Theorem 13.21 we can represent the 
matrix A + by (jU])- Therefore, we obtain for all i = 1, n, j = 1, s 



x 



i.) 



m E cdet i ((A*A). 4 (a* fc )) 

a Ik b kj = 2^ — 

k=l 



k=l 



E (A* A)? 



E Er=icdet,((A*A). J (a* fe ))^6 fcj 



E 



(A* A) 
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Since Yl a *k^kj 
k 



k 

Yl a 2k^kj 
k 



b j and denoting the j'th column of B by 



h j as well, then it follows (EH 



ii) If rank A = n, then by Corollary 13.11 A + = (A* A) 1 A*. Representing 
(A* A) -1 by we obtain for all % — 1, n, j = 1, s 



X j 



1 - 



ddetA 



k=l 



where is a left ijth cofactor of (A* A) for all i, j — 1, ...,n. From this by 
Lemma [2.21 and denoting the j'th column of B by bj, it follows (ITo"]) . ■ 

Corollary 4.1 (Theorem 3.1 in JTty ) Suppose 

AX = B (16) 

is a rig/it matrix equation, where {A,B} G M(n,H) are given, X G M(n,H) 
is unknown. If ddetA 7^ ; iaen (ESP /ias a unique solution, and the solution 
is 



cdet,(A*A) 



(17) 



Jl] ddetA 
where b j is the jth column o/B /or a// i,j = 1, ...,n. 

Definition 4.2 Consider a matrix equation 

XA = B, (18) 

where A G H mxn ,B G HP xn are given, X G M sxm is unknown. Suppose 

H L = {X| X G H sxm , ||XA - B|| = min}. 

Taen matrices X G HI sxm such that X G are called least squares solutions 
of the matrix equation ffify . If~X.LS = minxes ||X||, iaen X^g is called the 
minimum norm least squares solution of (Tl 



The following theorem can be obtained by analogy to Theorem 14.11 



10 



Theorem 4.3 The least squares solutions of (fig)) are 

X = BA+ + (I m - AA+)C, 

where C G HP XS is an arbitrary quaternion matrix and the minimum norm 
least squares solution is X^s = BA + . 

We denote BA* =: B = (%) G H sxm . 

Theorem 4.4 (i) //rank A = r < n < m, then for the minimum norm 
least squares solution X i>s = (xy) G H sxm of [Tfy) for all i — 1, ...,s, 
j = 1, m, we /iai>e 



E rdet^AA*)^.) 



E l(AA*) 2 | 



r. m 



(zij //rank A = m, then for X L5 = (xy) G H sxm o/ (EIP /or all i = 1, s, 
j = 1, m, we have 



rdet^AA*)^ (bj.) 
ddetA 



^; = ,~ r V 7 (2°) 



where bj zs t/je ith row o/B /or aZH = 1, s. 

Proof, (i) If rank A = r < n < m, then by Theorem 13.21 we can represent the 
matrix A + by (|T0|) . Therefore, for all % — 1, s, j = 1, m, we get 



E rdet,((AA*),>£ 
" h hk E l(AA-) S| 



a 

k.J I a 



aEl, 



r. m 



ELi hk E . rdet^AA*)^ 



Ct€lr,m{j} 



a 

k.J I a 



E l(AA*)sl 



a£l r 



Since E^fc a fc. = ( E ^fc a L E&ifc a fc2 "' E ^fe a L ) = bj and denoting 
the ith row of B by bj as well, then it follows ( TT9l) . 
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(ii) If rank A = m, then by Corollary 13.11 A + = A* (AA*) 1 . Representing 
(AA*) -1 by (jlj), we obtain for all % = 1, s, j = 1, m, 

1 - 

Xij = ddetA^^^' 

fc=i 

where i?^ is a left r/'th cofactor of (AA*) for all i,j = 1, ...,m. From this by 
Lemma [2. II and denoting the ith row of B by bj , it follows (|20|) . ■ 

Corollary 4.2 (Theorem 3.2 in [Tby ) Suppose 

XA = B (21) 

is a left matrix equation, where {A,B} G M(n, H) are given, X G M(n, H) 
zs unknown. 7/ddetA 7^ 0, t/ien ( fFZj) /ias a unique solution, and the solution 
%s 

rdet^AA*),, (b f .) 
%j ~ ddetA [ZZ) 

where bj is the ith column o/B for all i,j = 1, n. 

Definition 4.3 Consider a matrix equation 

AXB = D, (23) 

where A G H mxn ,B G H pxg ,D G H mxg are given, X G H nxp is unknown. 
Suppose 

H D = {X| X G H sxm , ||AXB - D|| = min}. 

T/ien matrices X G H sxm suc/i that X G ifo are called least squares solutions 
of the matrix equation l[23\) . //X L5 = minxes ||X|| ; then X L>S is called the 
minimum norm least squares solution of [23]) . 



The following important theorem is well-known. 

Theorem 4.5 (Tfflp The least squares solutions of l[23\) are 

X = A+DB+ + (I n - A+A)V + W(I P - BB+), 

where {V, W} C M nxp are arbitrary quaternion matrices and the least squares 
solution with minimum norm is Xls = A + DB + . 
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We denote D = A*DB*. 

Theorem 4.6 (i) //rank A = r\ < m and rankB = T2 < p, then for the 
minimum norm least square solution X^s = (xij) G M nxp of ( fl^j we 
have 

£ cdet,,((A*A). J (d«))J 

/3e J ri , 



or 



where 



• I j. 



£ (A*A)J £ l(BB*)°| 



E rdet^BB*)..^) 



:a*a)^ 



E l(BB* 



a6Z 



(24) 



(25) 



v «e-fr 2 ,p{j} 



rdet, (BB*). d fc . G HP , fc = 1, n (26) 



£ cdet.fl A A) , (d, 

. /3£j ril „{i} 



GH lxp , Z = l,...,p (27) 



are t/ie column vector and the row vector, respectively, d.;. and d j are 
the ith row and the jth column o/D for all % = 1, n, j = 1, 

fzzj //rank A = n one? rankB = p, then for X^ = (x^) G H nxp o/ 
we have 

cdet i (A*A). i (dj) 



.r 



or 



where 



x 



i.i 



ddetA ■ ddetB ' 

rdet^BB*),, (dfl 
ddetA • ddetB ' 



d* := rdet,(BB*),. d, 



,...,rdet,(BB*),, ( d n 



(28) 
(29) 

(30) 



d^ := (cdet i (A*A). < (da) ,...,cdeti(A*A). 4 (d. n )) (31) 
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are respectively the column-vector and the row-vector, is the ith 
row of D for all i = 1, ...,n, and d j is the jth column of D for all 
J = 1, -,P- 

(Hi) //rank A = n and rankB = r 2 < p, then for X^g = (zy) G HP xp of 
we have 



cdet,- (A* A) 



or 



ddetA £ l(BB*)°| ' 
E rdet^BB*) ,..((!£)) " 



(32) 



06^r 2 ,p{i} 



ddetA £ |(BB* 



(33) 



where d^ is 



an 



d d^ zs (E2P- 



f»J //rank A = r\ < n and rankB = p, then for X^s = G HP**' o/ 

rdet^BB*),. (dA) 



' .7 



E 

,n 



IA*A) 



• ddetB 



(34) 



or 



x 



E cdet^A'AJ.Jd*.)) 



J///3 



E 

w/iere d^. zs (C5J) and d^ zs ( flTj ). 



vA*A)^ 



■ ddetB 



(35) 



Proof, (i) If A G M™ xn , B G and r x < m, r 2 < p, then by Theorem 

321 the Moore-Penrose inverses A + = (ag) G e™ xm and B+ = (aj) G H 9Xp 
posses the following determinantal representations respectively, 

E cdet, ((A* A). . (a*.)) J 



, ^GJ ri ,„{i} 



E 



(A* A) 
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£ rdetj((BB*)j.(bJ)) J 



I* i a 

a I 



(36) 



E l(BB 

By Theorem 14.51 X^s = A + DB + and entries of X^s = (ajy) are 

q / m \ 



=1 \fc=i 

for all i = 1 , . . . , n, j = 1 , . . . , p. 

Denote by d. s the sth column of A*D =: D = (dy) G H nX9 for all 
s = 1, q. It follows from ^ a *k d ks — d.s that 



k=l 



E cdetj ((A*A) s (a* fc )) ^ 

/3€J ri ,n{i} 



k=l 



E (A* A) 

, 71 



E E cdet, ((A* A) . (a* fe )) J. 4, E cdet, ( (A* A) 2 ( d. 



E (A*A) 

, n 



E (A* A) J 



(38) 

Suppose e s and e. s are respectively the unit row- vector and the unit column- 
vector whose components are 0, except the sth components, which are 1. 
Substituting (1381) and (1361) in (|37|) . we obtain 



.r 



i.) 



/3eJ n ,n{i} 



E cdetj ( (A*A) ^ (d 



J E rdet,((BB*),,(b:.))- 

Q G/r 2 ,p{i} 



E (A*A)^ 



E l(BB*)Sl 



Since 



(39) 



z=i 



then we have 



x . 



15 



q p n 



EEE E cdet,((A*A). J (e. / ))^: t £ rdet,- ((BB*),. (e,)) I 



--it=ii=i/3e.J ri , n {i} 



Q Si"r 2 ,p{j} 



E (A* A) ^ E l(BB< 



p n 



EE E cdet i ((A*A).,(e.O)^« E rdet, ((BB*),.(e,)) 



i=H=l j g e j rin {i} 



«e/r 2 ,p{i} 



E (A* A) g E l(BB 



Denote by 



(40) 



4 



n 

cdet i ((A*A). i (d. t )) g = £ £ cdet,((A*A). 4 (e.,))^ 



( = 1 ,8eJ n ,n{i} 



the tth component of a row-vector = (d{^, d{^) for all t = l,...,p. 
Substituting it in (TJO]), we have 



Xij — 



E4 E rdet,((BB*),.(e,))- 

i=l ae-fr 2 , P {j} 



/SeJ, 



E (A*A)^ E l(BB*)« 



Q6/r 2 ,p 



Since ^ 4 e t = d^, then it follows 
If we denote by 



d U : =Z> E rdet,((BB*),,(e,)) 



£ rdet^CBB^.Cd,.)) « 



^e-fr 2 ,p{j} 



(41) 

the /th component of a column- vector d 1 * = (dfj, ...,df n ) T for all / = l,...,n 
and substitute it in (14131) . we obtain 



E E cdet J ((A*A). J (e. i ))^g 



— 



1=1 /3GJ ri ,„{i} 



j 71 



E (A* A) ^ E l(BB*)2| 



OS/r 2 ,p 
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Since £ e.jdg = d^-, then it follows (2 1). 



1=1 



A*A) _1 A" 



(ii) If rank A = n and rank B — p. then by Corolrary l3.1[ A + 
and B+ = B* (BB*) -1 . If we represent (A* A) 1 as the left inverse by ( EJ) 



and (BB*) 1 as the right inverse by (jlj) , then we obtain 



%21 %22 

\X n \ X n 2 

( d\\ d\2 

d<21 "22 



A*A i A*DB* BB H 



\d n i d, 



n2 



Xlp\ 

%2p 



X 



( T A J 1 



1 



ddctA 



np J 

4m 



dnm ) 



A 

'21 
tA tA 

Li 12 Li 



22 



1 



fR u R 

pB pB 

rt 12 ri 22 



ddctA 



2n 
B 

21 



\R^p R\ 



2p 



T A \ 

r A 

LnnJ 
rL p2 
R ^pJ 



where L A is a left ijth cofactor of (A*A) for all i, j = 1, 
right ijth cofactor of (BB*) for all i,j = 1, ...,p. This implies 

p 



n and Rfj is a 



X j 



E E L^dks 

s=l \k=l 



R 



(42) 



v ~ ddetA • ddetB ' 
for all i = 1, n, j = 1, ...,p. We obtain the sum in parentheses by Lemma 
and denote it, 



ks = cdeti(A*A) 



k=l 



-:d 



where d , is the sth column-vector of D for all s 



,p. Suppose dj 



(df\, . . . , d A p ) is the row- vector for all i = 1, n. Reducing the sum ^ d A R 



A pB 

•js 



by Lemma 12. 1\ we obtain an analog of Cramer's rule for the minimum norm 
least squares solution of (I2"3"j) by (I2"gj) . 

Interchanging the order of summation in (142 p . we have 



X ; 



k=l 



E ^ ( E d ks Rf 



v _ 



ddetA • ddetB 
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Further, we obtain the sum in parentheses by Lemma [2.11 and denote it, 

v 

J2dksRf s = rdet^BB*),. (d fc .) := , 

s=l 

where d^. is the kth row- vector of D for all k = 1, ...,n. Suppose d 1 * := 
(dfj, . . . , dnj) T is the column- vector for all j — 1, ...,p. Using Lemma l2~2l for 

n 

reducing the sum E L-j^dfj, we obtain an analog of Cramer's rule for f f23l) 
by d29]). 

(iii) If A G H™ xn , B G W* q and r x = n, r 2 < p, then by Remark O 
and Theorem 13.21 the Moore- Penrose inverses A + = (aj) G M axm and B + = 
(6^) G HF xp possess the following determinantal representations respectively, 

+ _ cdet,.(A*A), t (a*) 
°* ddetA 



K = „ D ^„, • (43) 



E rdet,((BB*),.(b*)) 
E l(BB*)2l 

Since by Theorem 14. 5[ X^s = A + DB + , then an entry of = is ([3j 
Denote by d. 5 the sth column of A*D =: D = G HP X9 for all s — 1, q. 
It follows from E a *fc4s = d. s that 



™ ™ r HetfA*Al (aM cdeti(A*A) . 

^ ^ ddetA ddetA 

fc=i fc=i 

Substituting (jUJ) and fj43l in (l3"Tj) . and using (|39|) we have 



9 cdet,(A*A 

E — 



a 



E rdet^BB'Ub:.)), 

a e/r 2 , P {j} 



ddetA E l( BB * 



E E E cdet,,(A*A) _ . (e.,)d,.6^ E rdet,((BB*V (e,)) 

s=lt=ll=l aG-fr 2 ,p{j} 

ddetA E l( BB *)al 
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E Ecdet,(A*A) .(e.,) d lt E rdet^BB*),-. (e,)) S 
ddetA E l(BB*)g| ' 
If we substitute (|4T|) in (1451) . then we get 

n 

Ecdet i (A*A). i (e.,) d, 



(45) 



n 

B 

_ 1=1 

Xij 



ddetA E l(BB*)2l ' 

n 

Since again 

E e./rfg = d B j} then it follows ([32]) , where d^ is (126]) . 
If we denote by 

d A ■= 

a it ■ 

n n 

^cdet,(A*A) A (d. t ) = J]cdet,(A*A) ..(e.,) 4 
i=l z=i 

the tth component of a row- vector d^ = (d A , d{£) for all i = l,...,p and 
substitute it in (j4"5j) . we obtain 

E4 E rdet,((BB*),.(e,))S 



ddetA E l( BB *y 



Since again E d^ e t. = then it follows ( 133|) . where d,^ is (l3Tj) . 

t=i 

(iiii) The proof is similar to the proof of (iii). ■ 

Corollary 4.3 (Theorem 3.3 in Tffl$ ) Suppose 

AXB = C (46) 

is a two-sided matrix equation, where {A,B,C} G M(n, Hf) are given, X G 
M(n,H) is unknown. //ddetA 7^ and ddetB 7^ , then (fffi| ) /ias a unique 
solution, and the solution is 

rdet,-(BB*),- (cf) , s 

JJ ddetA -ddetB ' 1 ; 
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or 

_ cdet t (A*A). t (c°) 
XlJ ~ ddetA-ddetB ' 1 8J 

where cf := (cdetj(A*A). 4 (c.i) , . . . , cdet i (A*A). i (c. n )) is £/ie row vector 
and := (rdet 3 -(BB*)j, (c 1 ) , . . . , rdet^BB*)^ (c n )) T is the column vector 

and Cj. ; c j are t/ie ii/i row vector and thejth column vector ofC, respectively, 
for all i,j = 1, n. 

Remark 4.1 In Eg. (24\), the index % in cdetj ((A* A) i (d^)) designates 

ith column of ((A* A) f (a*)), but in the submatrix ((A* A) . (a*-)) i the 
entries of a*, may be placed in a column with the another index. Similarly, 
we have for j in W3\) . In (1$ ) and / TTPj) we have equivalently. 



5 An example 



In this section, we give an example to illustrate our results. Let us consider 
the matrix equation 

AXB = D, (49) 



where 



A 



/ 1 


i 3 \ 


— k 


i 1 


k 


j -i 


u 


-1 i) 



B 



% 1 j 

j k -i 



D 



A i J'\ 

k i 

1 j 

\0 k i) 



Then we have 



/ 4 2z + 2j 2j + 2k 
A*A= -2i-2j 4 -2i-2k 
\-2j-2k 2i + 2k 4 

/2 + 2i + 2j -i + 2j - 2k\ 
D = A*DB* = 1 - i - j i-j-k . 

\ 1 - 2j 2% — k J 

Since ddet A = det A* A = and det A* A = (_ 2 ^_ 2j 2% \ =8^0 
then rank A = 2. Similarly rankB = 1. 
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So, we have the case (i) of Theorem 14. 61 We shall find the minimum norm 
least squares solution X iS of (l4Uj) by (I2~4"j) . We obtain 



E k bb * 

OfS/l, 2 



3 + 3 = 6, 



E (A* A) ^ 



det 

By (126|) . we can get 



/3eJ 2 . 3 
4 

2j - 2A; I 



2j + 2k 



det 
+ det 



2i + 2j 



+ 



'2 + 2z + 2/ 
1 - i - j 
l-2j 



2? - 2j 4 

4 -2? - 2A; 
2i + 2k 4 



-z + 2j + 2k^ 
i-j-k 
2% - k 



24. 



Since 



'2 + 2i + 2j 2i + 2j 2j + 2k 
(A* A).! (d^) = ( l-i-j 4 -2z-2k 

1 - 2j 2i + 2k 4 



then finally we obtain 



E odetftA^.^dB))! 



E (A*A)^ E l(BB* 



. /2 + 2z + 2j 2i + 2A , , /2 + 2« + 2j 2j + 2fc 
cdet, ( 1 _._. J 4 ^j+cdet^ 1-2j . ^ 

144 

A + 5i + Qj - k 
72 ' 



Similarly, 



Xl2 



. , -i + 2j -2k 2i + 2j\ , + /-i + 2j - 2A; 2j + 2k 

144 
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-1 - 2% + 5j - Ak 



72 



, 4 2 + 2z + 2 A . fl-i-j -2i - 2fc N 
cdet 2 n . . 1 . ) +cdeti ' 



^21 



^22 



2i-2j J ' 1 V 1 — 2 j 4 

144 

i - 2j - k 
72 ' 
2k^ 

-2i-2j i-j-k J 1 V 2z - k I 

144 

-2 + 2i + j - 



, 4 -z + 2j-2fc\ . fi-j-k -2% -2k 
cdet 2 . . . , + cdeti ' 



72 



4 2 + 2z + 2? \ . / 4 
cdet 2 . OI , n • + cdet 2 



£31 



^32 



2j -2k l-2j J z \2i + 2k 1 - 2j 

144 

l-4i- 3j 
72 ' 

-i + 2 j - 2&\ . x / 4 i-j-k 
C<1 ° t2 ( -2j-2fc 2i-j J +cdetl ^ + 2fc 2^-fc 

144 

3i - 3j - 2A; 



72 

Then 



1 



A + 5i + 6j-k -l-2i + 5j-Ak 
X L5 = ^[ i-2j-k -2 + 2i + j-k 
'- 1 l-4z-3j 3i-3j -2k 

is the minimum norm least squares solution of f l4"9]) . 
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